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ABSTRACT 

We present new integral-field spectroscopy in the outskirts of two nearby, edge-on, late-type galaxies to 
search for the Ha emission that is expected from the exposure of their hydrogen gas to the metagalactic ultravi- 
olet background (UVB). Despite the sensitivity of the VIRUS-P spectrograph on the McDonald 2.7m telescope 
to low surface brightness emission and the large field-of-view, we do not detect Ha to 5cr upper limits of 
6.4x 10-19 erg/s/cm^/D" in UGC 7321 and of 25x lO^i^ erg/s/cm^/D" in UGC 1281 in each of the hundreds 
of independent spatial elements (fibers). We fit gas distribution models from overlapping 21 cm data of HI, ex- 
trapolate one scale length beyond the HI data, and estimate predicted Ha surface brightness maps. We analyze 
three types of limits from the data with stacks formed from increasingly large spatial regions and compare to 
the model predictions: 1) single fibers, 2) convolution of the fiber grid with a Gaussian, circular kernel (10" 
full width half maximum), and 3) the coadded spectra from a few hundred fibers over the brightest model re- 
gions. None of these methods produce a significant detection (> 5a) with the most stringent constraints on 
the HI photoionization rate of r(z = 0) < 1.7 x 10"" s^^ in UGC 7321 and r(z = 0) < 14 x 10"" s^^ 
in UGC 1281. The UGC 7321 limit is below previous measurement limits and also below current theoretical 
models. Restricting the analysis to the fibers bound by the HI data leads to a comparable limit; the limit is 
T{z = 0) < 2.3 X 10^^'^ s~i in UGC 7321. We discuss how a low Lyman limit escape fraction in z ^ red- 
shift star forming galaxies might explain this lower than predicted UVB strength and the prospects of deeper 
data to make a direct detection. 

Subject headings: galaxies; evolution — diffuse radiation — intergalactic medium 



1. INTRODUCTION 

The strength of the metagalactic ultraviolet background 
(UVB) has great impact on theoretic al models of structure for- 
mation (e.g. lHaardt & Madaulll996l) and a variety of physical 
processes such a s the inhibition of small halo collapse (e.g. 
lEfstathioull 19921) . the intergalactic temperature and ionization 
state of the intergalactic medium (IGM) (e.g. iHui & Gnedir 
19971). and IGM metallicity determinations (e.g. lRauch et al 
1997al) . The likely contributors to the U VB are active galac 



tic nuclei and star formation in galax ies (ISchirber & BullocB 
120031; iFaucher-Giguere et alj l2008bl) with appear compati- 
ble w i th observed populations ( Gallego et al.lll995l ; iHopkinsI 
12004 iHopkins et al.ll200"7t iBouwens et al.1 120091) under rea- 
sonable corrections for dust attenuation, low luminosity ex- 
trapolations, redshift evolution, and escape f ractions. The 
strength of the UVB, especially at low redshift (iDave & TrippI 
1200 Ih is still highly uncertain despite its importance. Most 
recent efforts have focused on high redshifts, z > 2, 
where the strongest UVB measurements exist. For in s tance , 
the detailed history of star formation (iMadau et al.l 119991 ; 
iFaucher-Giguere et alJ l2008ah and the potential to measure 
individual active galactic nuclei (AGN) host h alo masses 
jLoeb & EisensteinI 19951: IFaucher-Giguere et al.ll2008cl) have 
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been explored. Measurements of the photoionization rate 
have used three methods: observations of Ha such as de- 
scrib ed in this paper, the line-of-sight proximity effect method 
(e.g. ICarswelletan[T98l iBaitUk et al.l figSSl). and the flu x 
decrement method (e.g. ICen et al.ll 1 994t iRauch e t"anil997bl) . 
The latter two require backlighting quasars and are therefore 
difficult or impossible at low redshift. We are motivated to 
constrain the current model with a different, low redshift mea- 
surement. Instead of using Lyman-a forest features, we pur- 
sue a measurement of the UVB powered. Ha emission that 
should occur in the outskirts of local disk galaxies. As a sec- 
ondary motivation, the kinematics of Ha at distances beyond 
HI data are importan t probes to the total da^rk ha lo masses in 
nearby disk galaxies ( IChristlein & Zaritskvll2008l) . 

Galactic disks are optically thick to Lyman limit 
photons and maintain their observed HI distributions 
through self-s hielding a gainst t he UVB. As recognized 
for decades (ISunvaevI Il969t iFelten & BergeronI 119691: 
iBochkarev & SunvaevllI977l) . the influence of the UVB may 
be investigated in the extreme outskirts of disks where the 
self-shielding begins to fail. These early works sought to mea- 
sure this effect through dis k truncation in HI . However, there 
appea r to be 
1993!) and 



cases with (Corbelli et al. 
without (Walsh et al. 1997; 



1989; vanGorko 
Carignan & Purto: 



1991 lOosterloo et al.ll2007l^ HI truncations above the critical 



column density predicted using current UVB estimates, im- 
plying that other processes may strip gas and mimic the result. 
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Moreover, reaching the UVB impUed truncation thresholds in 
21 cm measured HI would require rather long observations 
with current facilities. A more robust signature of the UVB 
strength would be the detection of the Ha in these outskirt 
regions. Ha has been found at such radii before in actively 
star fo rming and warped galaxies by iBland-Ha wthorn et alJ 
(119971) (hereafter BFQ) with Fabry-Perot staring measure- 
ments. However, the fi{Ha) = 2.3x10"^*^ erg/s/cm^/D" 
detection was interpreted to be due to non-UVB sources 
as indicated by an abnormally high [NII]A6548 to Ha ra- 
tio. Searches hav e also yielded li mits i n quiescent sy stems 
(IVogel etalJll995t IWevmann et aHIIOOll iMadsen et al.ll200l 
with an upper limit for the UVB photoionization rate, T, of 
r(z = 0) < 2.4 - 9.5 X 10"i-^s-i(2cr) being the deepest. 
The wide range due on this limit is due to gas cloud geometri- 
cal uncertainty. Despite the numerous theoretical implications 
and the efforts of numerous groups, a UVB powered Ha de- 
tection still awaits discovery. 

The tactical advantages we bring to this problem are deep 
surface brightness limits, a large two dimensional field of 
view through integral field spectroscopy compared to the pre- 
vious longslit and Fabry-Perot staring data, and target selec- 
tion of very high inclinations to maximize signal and min- 
imize contamination uncertainty. Our targets are edge-on, 
low surface brightness Sd galaxies that are rather isolated 
and minimally warped in order to avoid density distribution 
uncertainties and exposure to internally generated ionization 
from smaller radii. Indeeed, our most constraining target, 
UGC 7321, has a gas surface density below that required for 
significant star formation (Kennicutt.l989l) at all radii, as well 
as being unusually isolat ed with no known comp anions and 
minimal (< 3°) warping (lUson & Matthewsll2003b . 

In this paper we begin with a description of the simple ion- 
ization state and density model of disk galaxies that will be 
used to link a measured Ha surf ace b right ness with a partic- 
ular UVB photoiozation rate in In m.2\ we give disk 
param eter constraints based on fits to existing 21 cm data. In 
m.3\ we argue that UGC 7321 in particular is likely to ex- 
tend its HI profile beyond the current 21 cm limits without 
truncation. In addition, the HI observations of UGC 7321 are 
amongst the most sensitive such measurements published to- 
date. The 21 cm data allow a very precise model to be made 
for the gas distribution in the galaxy outskirts at the locations 
where we search for Ha emission. Next, in S|3] we present 
deep integral field spectroscopy observations at radii corre- 
sponding to the outermost detections of 21 cm emission and 
beyond. We describe the choices made to stack spectra on 
various spatial scales. The stacked spectra are searched for 
Ha detections and upper limits are derived. Particular focus 
is given to systematic errors. Finally, in 21 we discuss the 
context, the likely cause of the unexpectedly low limit, and 
further observations that can confirm our conclusions. The 
Appendix lAl provides the analytic details necessary to con- 
struct the full and general Ha surface brightness distribution 
model. We will quote most of the surface brightness limits 
in units of erg s^^ cm~^ arcsec^^, but for easy comparison 
to alternative units we note the conversion at the wavelength 
of Ha of 1 millirayleigh (mR)= 5.66 x 10~^^ erg s^^ cm~^ 
arcsec"^ = 2.8 x 10~^ cm~^ pc in emission measure assum- 
ing the case B coefficient we adopt. 

2. HI BASED MODELS AND Ho PREDICTIONS 
2.1. Model assumptions 



A three dimensional gas density distribution must be in- 
ferred in order to translate Ha surface brightness into a 
UVB strength. BFQ made estimates assuming exponen- 
tial forms both radially and vertically in the gas distribu- 
tion with a plane parallel assumption. Motiv ated by the reg- 
ular HI stru cture on local scales (Garcia-Ruiz et al.l I2002t 
lUson & Mat thews 2003) of our chosen targets showing sim- 
ple exponential trends and needing an extrapolated model in 
gas density for interpretation of UVB limits, we also assume 
exponential forms. 

In order to interpret Ha measurements generically inside 
and outside of the UVB photoionization front around gaseous 
disks, we have generalized the model of BFQ. Some toy cal- 
culations in the model also show the importance of high in- 
clination selection to make the deepest possible UVB con- 
straints. This high inclination boon has been known be- 
fore, but not carefully followed in earlier works' target se- 
lection. The model assumes both regular gas distributions 
and sharp photoionization transitions in a plane parallel ap- 
proximation under arbitrary disk inclinations and sight lines. 
Our model assumes sharp photoionization fronts exist. We 
verify this assumption by estimating the Lyman limit pho- 
ton mean free path at the midplane ionization front. In their 
Equation 3, BFQ estimate the hydrogen density at this point 
as uh ~ 0.05cm^'^. The Lyman limit photon mean free 
path is given by Imvf ^ ( n x Oi/)^^ ~ l.lpc with Qi, 
jOsterbrock & Ferlandl 120061) as the hydrogen Lyman limit 
photoionization cross section. This is much smaller than 
the common disk scale lengths in either direction. The ver- 
tical scales for cold disk galaxies are of order 1 00 pc or 
greater More sophisticate d models can be made (iMalonevI 
119931; iDove & Shulilll994l) by solving for the ionization and 
excitation states of hydrogen and helium with full radiative 
transfer solutions in a grid of plane-parallel gas layers, but 
such an analysis is beyond the scope of this work. 

The forthcoming derivation follows BFQ equations 1-6. 
The important differences are that this derivation is general- 
ized for any viewing inclination, i, and for arbitrary position- 
ing of the spectral data in the galaxy's field of observation. 
The BFQ derivations were specifically for i = 0° and the field 
position along the major axis where all gas is photoionized. 
We denote the generic surface brightness in Ha as ji. We de- 
note /xg as the special case of the peak Ha surface brightness 
where the photoionization front intersects the disk midplane. 
Our results reduce to the BFQ values of i^lq for i = 0°. In 
Equation [U we give the assumed gas distribution in cylindri- 
cal coordinates R and z with radial scale-length h^, vertical 
scale length h,, and central hydrogen density uq. 

z) = noexp(-^)exp(-^) (1) 

The commonly assumed form of the UVB spectrum is given 
in Equation|2]where ly is the frequency, vq is the Lyman limit 
frequency, Jq is the UVB strength at the Lyman limit in units 
of erg cm^^ s^^ Hz^^ sr^^, and f3 is the UVB spectral index. 

j.^joC^y (2) 

Another common form of quoting the UVB strength is with 
the UVB photoionization rate, F. We show this form in Equa- 
tion[3] where h is Planck's constant, (t(i/) is the hydrogen pho- 
toionization cross section, and a^=(j{i'o) is the Lyman limit 
cross section. The final equality in Equation[3]comes from the 
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standard power law approxim ation to the cross section shape We use ct 



stanaara power law approxm 
jOsterbrock & Ferlan(3i200"6l) 



r = 47r 



hv 



dv ~ 



hx{-i + l3) 



(3) 



In Equation |4] we equate recombination and ionization rates 
under a plane parallel approximation. For the radial regions 
where any self-shielding can take place, we consider the top 
and bottom of the disk to each see incident flux from only half 
their total solid angle. We define rig as the electron density. 
Up as the proton density, ^ as the ionization fraction, as as 
the case B recombination coefficient, and Zc{R) as the height 
above the midplane to which the photoionization front pen- 
etrates at radius R. We define e as the volume filling factor, 
assumed to be spatially invariant. A clumpy gas distribution 
can, to first order, be represented by using this term somewhat 
lower than the nominal value of unity. With the assumption of 
sharp ionization boundaries, we can equate the gas densities 
as He = Up — S^uh at radii beyond the photoionization front. 



/>oo I^C 

as / £_'^enj^{R, z)dz ~ I 



2'KJ,y 27rJo 

-dv = . „ (4) 



hv 



h/3 



We next define a threshold radius, Vc, to which the UVB pen- 
etrates fully through the disk plane, so zdvc) = 0. The solu- 
tion of Equation|4]leads to Equations |5] and |6] 

r, = i\n{2^haBa,nlh,P) - ln(r x (3 + /?))) x K/2 (5) 
±(r, -R)x hJK ■■R<rc 



Zc{R) 



0: R>rc 



(6) 



Next, we define the variable p as the distance from the disk's 
midplane along the line-of-sight, spanning —oo to the ob- 
server and oo away from the observer We also define the 
major axis position bi, and minor axis position 62 as the ob- 
served field positions projected onto the sky. Finally, we rep- 
resent the galaxy's inclination with i. Simple transformations 
to cylindrical coordinates give the expressions in Equations |7] 
andg] 



2p&2 sin i cos i (7) 



R = 



sin^ i 



^2 cos^ i + 6^ — 2pfo2 sini cosz (8) 



The Ha surface brightness, /i(6i,62), follows directly from 
a line-of-sight integration. The full evaluation of /i(6i,fe2) 
involves finding the values of p that intersect the photoion- 
ization surface described by Zc{R) with field positions bi and 
62- The analytic solutions to those intersections are given in 
Appendix |A] The solution for the special case at field po- 
sition 62 = and 61 = Tc gives the aforementioned peak 
surface brightness pQ, which is itself a useful measurement 
parameter, as related in Equation |9] In Equation |9] 7 is the 
generally non-analytic integration of the emissivity along the 

line of sight, a|J^ is the case B effective Ha recombination 
coefficient, and il is the full sky solid angle of 4tt sr For 
t he g as density parameters we derive in our target galaxies in 
j i2.2l and the areas we observe in fj3] the face-on column den- 
sities of total hydrogen are sufficient (> 10^^ cm~^) to be 
everywhere optically thick to Lyman limit photons, let alone 
Lyman-/? and the other important lower energy transitions. 



eff 

He 



1.17 X 10 cm^ s ^ as appropriate for 



T= lO^'K (lOsterbrock & Ferlandll2006h . 



2hi^HaQ:^Ji /o°° ne{R, z)np{R, z)dp 



fJ-Q 



2(^ea^^l^nlhvHal 
with 7 = 



, 2pcosi 2Wr2 +p2sin i 
exp( r 7 )dp 

'^2 "-r 



(9) 



We explain the use of certain constants and assumed values 
to Equation |9] The ionization fraction is assumed to be unity 
by the earlier discussion of the Lyman limit mean free path. 
The volume filling factor may approach unity as there is no in- 
dication of star formation at extended scales in these galaxies. 
We will discuss the evidence for the absence of extended star 
formation in Sect ion IH Furthermore , the d eprojection of the 
HI distribution in lUson & Matthewsl (l2003h gives a peak sur- 
face density of only 5.8 M© pc"^ at the center of UGC 7321. 
The surface density drops by over an order of magnitude at the 
locations we observe. These surface densities are well below 
the dy namical criterion for efficient star formation (iKennicutll 
and make a smooth gas distribution plausible. It is not 
possible to exclude small scale dumpiness, so we retain the 
volume filling factor The case B and Ha effective recom- 
bination c oefficients are d e pende nt on electron temperature. 
Following I Wevmann et al.l ( 1200 lb and the discussion therein, 
we adopt T= 10,000 K and the values oflOsterbrock & Ferlandl 
(^06). The true electron temperature may plausibly be dif- 
ferent by a factor of two, leading to corresponding changes in 

a|/^ and as of the same order of magnitude. However, the 
linearization in F of Equation |9] makes the surface brightness 
depend on the ratio of these two recombination coefficients, 
so their similar behavior with electron temperature cancels. 
For consistency with previous works, we do not propagate the 
recombination coefficient uncertainties as systematics to the 
final UVB limit. 

Some brief numerical examples set the expected orders 
of magnitude, quantify the achievable limits under different 
galaxy geometries, and illustrate the important parameter de- 
pendencies under linear expansions. We look at some trial 
cases with e = 1, K ^ 1000 pc, = 100 pc, /3 = 1.8, 
r = 4 X 10"" s-\ and no = 5 cm^^. For i = 0°, 
7 = cxp(-2rc//ir) X h^/2 so /x « 3.0 x lO'^" erg/s/cm^/D". 
For i = 90°, 7 = rcKi{2rc/hr). Ki{x) is the modified 
Bessel function. In this case, p ss 1.3 x 10"^* erg/s/cm^/D". 
For this work's applications, the surface brightness profiles 
are smoothed by seeing and sampled by large fibers. Re- 
alistic smoothing and sampling, of order several arseconds, 
can lower these peak values by several tens of percent. We 
will assume for all c alculations that (3 = 1.8 as motivated 
by previous models dShull et al.lll999t) and to aid the com- 
parison with previous observationa l work that used the same 
assumption (I Wevmann et al.ll2001h . We note that pQ scales 
exactly linearly with F when viewed face-on and nearly lin- 
early for all other inclinations. This is easily demonstrated 
by taking the large argument asymptotic behavior of the mod- 
ified Bessel function which yields a linear scaling in F af- 
ter a first order expansion. We show the small error caused 
by assuming a linear relation between ^0 and F in Figure 
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[U for reasonable geometries. All further estimations of F 
in this work will be made in the linear approximation. We 
have linearized our estimate around F = 4 x 10~^^ s~^ 
because we consider it the bes t curre nt estimation from the 
work of iFaucher-Giguere et al.l (l2009h . However, any initial 
value would have worked as the only effect of a particular 
choice is that the small non-linearities pivot around the sim- 
ulation UVB choice, but this error is negligible compared to 
our other error terms. The discussed numerical example be- 
tween i = 0° and i = 90° also shows how the selection 
of thin, edge-on disks can exploit a particular flux limit to 
a (30 — 50x) stronger UVB constraint than for face-on disks. 
We also emphasize with Equation|9]that the first order effects 
near z=90° on distance, volume filling factor, and gas density 
cancel out; /ip only has first order dependence on i, the ratio 
of scale lengths, and F. 

2.2. HI data 

Three-dimensional gas distributions must be inferred for in- 
dividual galaxies to interpret Ha surface brightness and to 
guide the stacking choices amongst fibers. We will use such 
fits to extrapolate the density profiles to larger radii where 
the Ha emission is predicted to reach peak surface bright- 
ness. The parameters from stellar distributions could po- 
tentially be used, but 21 cm measured HI is the more rele- 
vant indicator to ioni zed hydrogen. We adopt distances of 10 
Mpc for UG C 7321 (lUson & Matthewsl l200l and 5 Mpc for 
UGC 1281 (iTuUy et al.ll2006h . Low redshift surface bright- 
ness is insensitive to distance, so the exact distances are unim- 
portant to this work. Different literature estimations have 50% 
and 10% rms ranges for the UGC 7321 and UGC 1281 dis- 
tances respectively. We indicate scale lengths by the terms 
dio as the actual distance to UGC 7321 in units of 10 Mpc 
and as the actual distance to UGC 1281 in units of 5 Mpc. 
For reference, the scale conversions become 48.5(iio pc/" for 
UGC 7321 and 24.2^5 pc/" for UGC 1281. 

UGC 7321 was observed by one of us in collaboration with 
L. D. Matthews (JJson & Matthews 20Q3|) using the second 
most-compact (C) configuration of the VLA^ which includes 
some of the shortest spacings available and their full cover- 
age, deep observations yielded spacings down to 28m, close 
to the dish diameter Their quasi-naturally weighted ("ro- 
bust" parameter TZ = +1) images recovered the full single- 
dish flux and, moreover, their single-dish equivalent spec- 
trum matched the features of the best single-dish spectrum 
to wit hin the (higher) un c ertain ty of the single-dish observa- 
tions (lUson & Matthewsl ( 120031) . fig. 6). For this paper, we 
have used their quasi-uniformly {TZ = —1) weighted images 
because of their better resolution (~^12" ) although the some- 
what higher (45%) noise level only recovers ^96% of the total 
flux. However, the five parameter model fits to the zeroth mo- 
ment maps, described below, recover some of the lost flux and 
the remaining uncertainties only slightly shift the position of 
the predicted Ha peak. 

For UGC 1281, we have reduced the raw data from the 
VLA archive. It was observed under proposal AZ097 on 
1997 December 26 in the most compact (D) configuration 
for a total of 3 hours on source with interspersed observa- 
tions of the strong, primary calibrator J0137H-3309 (3C48) for 
which we have adopted the VLA recommended flux density 

^ The Very Lai'ge Array of the National Radio Astronomy Observatory 
is a facility of the National Science Foundation, operated under cooperative 
agreement by Associated Universities, Inc. 



of 15.9 Jy. The observations were spaced over a range of 
±3 hours in H.A. giving excellent uv-coverage and images 
with 127 channels of width ^2.6 km s after standard on- 
line Hanning-smoothing. The angular scale that corresponds 
to the shortest baseline (~900") is sufficiently larger than the 
largest single-channel galaxy extent (~285") that the array 
should have recovered the total HI flux. We followed the 
same reduction steps as for UGC 7321 (lUson & Matthewsl 
12003b to obtain a "cube" of spectral images using nearly- 
natural weighting (TZ = +1) which gave images with reso- 
lution ~'51" which were free of artifacts to the rms sensitivity 
(T ^1.0 mJy/beam per channel. We computed moment maps 
after applying a standard "1-cr cutoff" evaluated on a cube 
Gaussian- smoothed spatially to 70" and Hanning-smoothed 
in frequency which led to a total HI flux of 41 ± 2 Jy km s^^ 
corresponding to a mass of 2.3 x 10'* d^MQ. The total flux 
is in good agreement with the values in the literature which 
range from (35.8 to 38.9) Jy km/s from two different single- 
dish measurements (iHuchtmeierl 1 1 9891) with the spread and 
uncertainty due in part to some ringing from strong in-band 
HI emission from the Milky Way as well as to calibration un- 
certainties. Again, we have obtained a spectral "cube" with 
nearly-uniform weighting {TZ = —1) which gave images with 
resolution of ~'42" with rms sensitivity cr ~ 1 .5 mJy/beam per 
channel. As in UGC 7321, the higher noise level results in a 
slightly lower total flux, 39 ± 2 Jy km s~^. 

Next, we characterize the HI distributions nearest our 
Ha observations. We have derived five parameter fits in 
no,hr, h,, i, and position angle to the zeroth moment maps 
of UGC 7321 and UGC 1281 through non-linear least squares 
minimization. The models include convolution to the instru- 
mental beams of ~ 1 2" and ~^42" FWHM and sampling of Ha 
appropriate to the fiber data. Bot h the maps show at lea st two 
major axis power law slopes, as iChris tlein et al.l (120101) have 
found to be common in extended gaseous disk gas. We do 
not try to model the full gas distributions, but only the large 
radius trends by restricting the fits heavily to the outermost 
data regions. Still, the model fits deviate from the data by 
an amount that exceeds the observational errors. Some minor 
warps and substructure are visible. The formal errors in the 
total line intensity images are 15 x 10^^ cm^^ for UGC 7321 
and 5 x 10^® cm^^ for UGC 1281, which are both far smaller 
than the residuals to the best fit models. In order to capture 
the systematic model errors, we have made Monte Carlo sim- 
ulations between the data and the best fit models to create 
68% confidence intervals as given in Table [1] for all disk pa- 
rameters and Ha observables. The perturbations in the Monte 
Carlo simulations are made from the residuals of the best fit 
model, not the statistical errors, to include the influence of 
systematics. These simulations allow us to create three types 
of Ha surface brightness prediction, with different scales of 
spatial co-addition, under an assumed F. Note that many 
of the individual disk parameters in Table [1] have large rela- 
tive uncertainties, but the surface brightness predictions have 
small relative uncertainties. The disk parameters share de- 
generacies, as captured in the Monte Carlo simulations, that 
create highly certain Ha predictions despite the individually 
uncertain gas parameters. Predictions can be made for indi- 
vidual fibers, but to both mitigate the model uncertainties and 
improve our limits, we include predictions with a 10"xlO" 
FWHM convolution sampled near the peak surface brightness 
positions. The exact choice of kernel size is not important, 
but is chosen to combine several neighboring fibers. Finally, 
we include a prediction for the average surface brightness of 



UVXHa 



5 



all fibers expected to sample ji > 10~^^ erg/s/cm^/D". These 
various predictions will be compared to co-added data in |J3] 
We give in Figure |2] the HI fits along major and minor axis 
cuts. The fits to UGC 7321 use all the HI data beyond an in- 
ner radius cut, which was chosen to avoid a substructure bump 
near R w 140". The fits to UGC 1281 are more constrained 
with both an inner and outer radius cut. The outer cu t is to 
exclude a known ~ 8° warp dOarcfa-Ruiz et al.l 120021) . The 
fitting function assumes a single position angle at all radii and 
does not describe warps well. We have investigated the disk's 
outer behavior by also deriving fits from the R>220" data 
alone. A radial scale length compatible with, but noisier than, 
the Table [T] value was found with a significant change in posi- 
tion angle, a reflection of the warp. 

2.3. HI bounded limit 

While we cannot definitively prove that these galaxies 
maintain their extrapolated hydrogen profiles over the galac- 
tocentric distances we will discuss in i)2.2| such an assump- 
tion is reasonable with the current 21 cm data. There is no evi- 
dence for flaring in these galaxies, and the s uperthin shape im- 
plies a n undisturbed history. In UGC 7321. IUson & Matthewsl 
(120031) have searched for low-mass compansions and found 
none to the limit of Mhi = 2.2 x 10^ Mq within 12'(35 
kpc). The nearest optical companions are two dwarf galaxies 
at minimum distances of 340 kpc, implying minimum times 
to last encounter of 1.6 x 10^ years. So, it is unlikely that gas 
has been stripped from the regions over which we have ex- 
trapolated a density profile. However, we have calculated an 
alternative limit using data bounded by the HI d a ta in a man- 
ner similar to the analy s is in[S tpcke et all (1 1 99 Ih : iVogel et al.l 
([T995h : iDonahue et all (119951) ? I Wevmann et al.l (12001 " as an 
alternative, which is equivalent to assuming that the gas is 
completely truncated where the 21 cm signal falls below the 
noise. In those works, a single, simple equation based on 
global photoionization equilibrium is used and here repeated 
in Equation [To] 

^ ^ j,3 + /3 ^ iHa Aproj ^^^^ 
4ai.^ fafna Aot 

The variable $ is the one-sided incident ionizing UVB flux 
in units of cm^^ s^^, Ihu is the Ha surface brightness in 
units of /iR, fa is the fraction of incident photons that become 
absorbed when passing through the face-on cloud, fHa is the 
fraction of excited recombinations that produce an Ha pho- 
ton, Aproj is the projected area covered by spectroscopy and 
21cm data, and Atot is the total surface area for the regions 
in projection that can absorb Lyman limit photons. The area 
aspect ratio is usually determined from 21cm data. This cal- 
culation takes no account of the spatial stratification between 
21cm and Ha that can realistically occur for very thin gas dis- 
tributions, as we will see later in ij2.2l where the predicted Ha 
surface brightness is derived, and requires Ha searches and 
interpretations to be restricted to area covered by deep 21cm 
data. However, for mild aspect ratios (^< 10) or large 21cm 
beams, this method delivers similar predictions as those in 

We now discuss the evaluation of the few terms in this 
model. The assumption in the HI bounded limit is that the 
hydrogen resides within some well-defined area represented 
by the noise floor of the 21cm data. It is not obvious how the 
area should be defined in a continuous gas distribution, but we 
adopt the photoionization front we have previously defined in 



Equations |5] and |6] as a realistic edge. In ^2.21 we determine 
gas geometries for our target galaxies. In particular for the 
area in UGC 7321 covered by fibers, with Nhi > 10^^ cm"^, 
and the parameters in Tabled we find (xTT") '^'^■^-i.t 
This value is in good agreement with the 21cm axis ratio of 
29 det ermined at the 10^" cm^^ contour in lUson & Matthewsl 
( I2003L Table 3). By adopting this distribution in face-on col- 
umn density and a UVB spectral index of /3 = 1.8, we can 
evaluate fa- We find {^^) = 22.8j:i:8 UGC 7321. 
With the same calculations applied to UGC 1281, we find 
{ ±±) = ^Q-O- ii and (^) = 13.61^2:^. Identically 
to lWevmann et al.l ( 1200 Ih . we adopt fna = 0.45 as appropri- 
ate for case B and a lO^K electron temperature. We also carry 
out this analysis in Tables [T]and|2]for continuity with previous 
work, but we emphasize that our preferred limit comes from 
the comparisons to the model in as it incorporates the 
spatial segregation between the brightest Ha regions and the 
HI data that is natural in very thin, edge-on geometries. 

3. Ha DATA AND ANALYSIS 

We have obtained new integral field spectroscopy posi- 
tioned along the major axes of UGC 7321 and UGC 1281 
targeting Ha with the Visible Integ ral-field RepUca ble Unit 
Spectrograph Prototype (VIRUS -P, iHill et al.l l2008h on the 
McDonald 2.7m telescope. We observed UGC 1281 on 2009 
October 22-24 with R = 1288 from 4700-6990A for 21 pho- 
tometric hours and UGC 7321 on 2010 April 9 and 11 with 
a resolution of R = 3860 from 6040-6740A for 15 hours 
under non-photometric conditions. Between the R = 1288 
and R = 3860 observations, made possible by a new grating, 
we not only gain in sensitivity scaled by the square root of 
the resolution but resolve the bright skylines, OH A6568.779 
and geocoronal Ha, from our target wavelengths. We have 
set the controller to bin pixels by two in the wavelength di- 
rection which samples the spectra just at the Nyquist cri- 
terion and minimizes read noise. The VIRUS-P field cov- 
ers a I '.6x1 '.6 field with 246 fibers of 2'.'05 radius with a 
one-third fill-factor. We split our observations into three 
dithers to cover the entire field. In UGC 1281 we split 
our time further between two overlapping fields to cover the 
outer plane better in the presence of a possible <8° warp 
{Garci'a-Ruiz et al. 200^ yielding a tota l of six dithers . Spec - 
trophotometric flux standard stars from iMassev et alj ( 119881) 
were measured once or twice nightly. We tracked the trans- 
parency through the offset guiding cam e ra. Galactic ex- 
tinction corrections (ISchlegelet al.1 119981: IO'Donnelll[T99l 
were made with Av=0.09 and Av=0.15 for UGC 7321 and 
UGC 1281 respectively. A spectral airmass/extinction curve 
specifically modelled for the McDonald Observatory site was 
applied. We estimate its systematic uncertainty by compar- 
ing it to the Kitt Peak curve supplied with the IRAF pack- 
age onedspec. We find a 20% rms difference between the 
wavelengths of 6000-7000A. The two curves deviate system- 
atically at A > 5900A. We believe the site specific McDonald 
curve to be more accurate to our data. However, we propagate 
the difference as a potential, systematic uncertainty. The flux 
calibration uncertainty due to the airmass/extinction curve at 
the data's median airmass of 1.09 is ±0.023 magnitudes. 

3.1. Flux calibration 
The 8' offset guiding camera is an Apogee Alta with a 
20.25n' field-of-view under a Bh-V (Amean = 5000A) filter. 
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Guider images were read out and saved every few seconds. 
Stacks of guider images that overlapped in time with each in- 
dividual VIRUS-P exposure (20 minutes each on UGC 1281, 
30 minutes each on UGC 7321, and 1 minute each on the flux 
standards) were combined. We make a relative photometry 
correction to each science frame based on the stack of guider 
images taken simultaneously with the VIRUS-P data. Typi- 
cally, ten stars per guider frame were available for photome- 
try. 

We have switched from the standard stars to the science 
targets with gaps of less than 5 minutes and assumed the con- 
ditions to be constant over that time and between the standard 
star and galaxy positions to make the absolute flux calibration. 
The observations of standard stars were taken during the most 
photometrically stable periods during each night to mitigate 
this potential source of error. Even so, the final flux calibra- 
tion factor we apply may have systematic errors. We assess 
this error by considering the 5 observations of 2 standards, 
PG1708H-602 and Feige 34, taken along with the UGC 7321 
data and the 3 observations of 1 standard, Feige 110, taken 
along with the UGC 1281 data. The distribution in flux cal- 
ibrations is wavelength-independent over our observed range 
with a 6.8% rms and 2.2% rms respectively. These estimates 
also capture possible variation in transparency with on-sky 
position. They are reported in Table [2] along with the pos- 
sible error in the extinction curve between the effective wave- 
length of the guider and the wavelength of Ha. For the non- 
photometric data on UGC 7321, we measured a median zero- 
point change, Azp, of 0.276 magnitudes and a 68% range of 
0. 171-0.382 magnitudes over the two nights. The more nearly 
photometric data on UGC 1281 had median Azp =0.057 
magnitudes and a 68% range of 0.043-0.077 magnitudes over 
the three nights. 

3.2. Sky background subtraction 

The choice of sky subtraction is particularly important for 
this work which reaches for flux limits far below the average 
sky brightness. If the science field were covered with source 
emission, sky nods would be necessary. Then, the time vari- 
ability of the OH and geocoronal Ha sky lines would form 
important systematic error sources. Fortunately, the large 
VIRUS-P field of view and selection of extremely thin, edge- 
on target galaxies affords a subset of fibers that contain a neg- 
ligible amount of source flux to serve as simultaneously mea- 
sured sky fibers. We selected fibers sufficiently far from the 
major axis such that the models predicted /i < 2 x 10~^^ 
erg/s/cm2/n"(with the basehne T = 4 x lO^i"^ s^^), or lOOx 
below the expected peak surface brightness, to be used for sky 
subtraction. This cut left 24% and 44% of the fibers for sky 
estimation in UGC 7321 and UGC 1281 respectively. We ex- 
perimented with moving this sky fiber cut up and down by a 
factor of five and found no difference in the final upper lim- 
its to the UVB strength. Depending on the number of fibers 
co-added, the statistical Ha flux errors presented here reach 
to 300 X dimmer than the sky level. Without simultaneously 
measured sky background, the systematics of sky nods would 
quickly dominate the limits. 

3.3. Data reduction 

The data reduction, optimal background subtraction, and 
search for emission lines were complete d with algorithms de - 
veloped for a Lyman-a emitter survey ( lAdams et al.ll2010l) . 
We summarize here the important steps. First, overscans and a 



master bias frame are subtracted from each frame. The wave- 
length solution for each fiber is fit as a fourth order poly- 
nomial to ^ 30 emission lines from HgCd lamps passing 
through the entire telescope light path. The residuals to the 
solution are of order one hundredth of a resolution element. 
Flat fields precise to < 1% are made from tw ilight flats with 
the solar spectrum removed by a b-spline fit (lDierckxlll993h 
and division. This fit method is the same as we apply to fitting 
and subtracting the sky background and has important advan- 
tages over data interpolation. By avoiding data resampling, 
we keep the errors largely uncorrelated. Small distortions of 
the instrument camera over a regular pixel grid lead to the 
spectrum from each fiber being sampled at slightly different 
wavelengths. By considering a collection of fibers together in 
a fit, the spectrum is oversampled, and we can recover nearly 
blended features. This method delivers an optimal spectral 
model robust against cosmic rays and without the residuals 
that linear interpolation can create. A thorough description 
of b-spline fits as applied to astronomy datasets can be found 
in iKelsoiil (l2003h . The next step in the data reduction is to 
fit and subtract a b-spline sky background modelled from se- 
lected sky fibers. Next, cosmic rays are masked by finding 
all pixels that deviate from the other pixels in the same fiber 
by some large threshold value. Some dim cosmic rays are 
missed by this step, but are rejected when combining multiple 
frames. We have chosen a threshold that misses the weakest 
^•20% of cosmic rays for direct masking in this work. The 
exact threshold does not affect the results. The frame is then 
flux calibrated with the non-photometric zeropoint correction 
and airmass correction applied. Finally, a one dimensional fi- 
nal spectrum for each fiber position is created by combining 
all the frames taken at the same dither position and running 
across the 5 pixel cross-dispersion aperture. For the final es- 
timate to be immune to remaining cosmic rays we have used 
the biweight estimator (Beers et al. 1990) at this step. Our 
pipeline makes no cross-talk correction since we restrict our 
cross-dispersion apertures to 5 pixels where the fiber separa- 
tions are typically 8 pixels and the cross-dispersion FWHMs 
are typically 4 pixels. This leads to, at most, 10% contami- 
nation from neighboring fibers and becomes especially trivial 
when considering large collections of fibers as an aperture. 
The scattered li ght propertie s of th e instrument have been 
characterized in lAdams et al.l (120081) and, particularly at Ha 
wavelengths, no scattered light or ghost patterns are found. 
The spectral resolution varies by < 5% for all fibers at a com- 
mon wavelength due to careful design and alignment of the 
spectrograph camera. We have made no corrections by con- 
volution to a common resolution. The effect of the spectral 
resolution variation and the background subtraction scheme is 
to leave residuals under bright s kylin es. We characterize the 
spectral resolution systematic in 8 13.51 Given the large number 
of independent spectral elements in VIRUS-P data (126,000 
in each dither), we must choose a high significance cut. At 
5(7 significance, the chance of noise leading to a detection at a 
particular wavelength in a particular dither is only 1 in 14,000. 
We choose to quote this limit as sufficiently conservative. 

3.4. Emission line detection 

We describe here an automated emission line search algo- 
rithm to work with a sky background and continuum sub- 
tracted spectrum or stacks of spectra. By applying this search, 
we robustly find all significant emission lines at all redshifts. 
In practice, we find no significant Ha emission with plausi- 
ble velocity offsets in any fiber for either galaxy. Plausible 
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velocity offsets are determined by the HI rotation curves. In 
UGC 7321, for example, the rotation curve is flat over our data 
range with variations of only ±10 km s~^. The gas disper- 
sion is measured in the HI data to be near 7 km s~^ subject to 
the li mitation of the 5 km s^^ resolution jUson & Matthewsl 
l2003h . Over a very conservative ±100 km s~^ (2. 2 A) range 
around our target wavelengths, the flux limit is flat. First, 
spectral pixels at any wavelength that exceed the noise by Ic 
are treated as seeds. Around each seed, we fit Gaussian pro- 
files of variable intensity, width, and central wavelength. The 
S/N of an emission line is then calculated by summing all pix- 
els and errors in quadrature within ±2(Tres for the wavelength 
of interest where ares is the instrumental dispersion. In the 
UGC 7321 data, ares =33 km s^^, and in the UGC 1281 data, 
ares =100 km s^^. When quoting limits on undetected emis- 
sion lines, we sum in quadrature the errors within the same 
spectral window. These steps in error combination consider 
both the statistical errors in the reduced data and the system- 
atic error based on ill-m atched spectral resolution between 
fibers discussed in ^13.61 A spectral correction factor is di- 
vided into the detections and limits to consider the fraction of 
a Gaussian function's flux that fall s outside of the considered 
window as f^pec = eif {V^ares/ \/cr^es + ^det) where adet is 
the detected emission line width. This same factor determines 
the degradation in flux limit for broad line detections. For un- 
resolved limits, adet is considered to be zero and the spectral 
correction (ffp\c) evaluates as xl.05. In practice, we make 
no significant detections within ±500 km/s of the HI based 
expected velocity in either galaxy. The average HI heliocen- 
tric velocities of UGC 7 321 and UGC 1 281 are 407 km/s 
(lUson & Matthewsl l2003h and 157 km/s (iGarcfa-Ruiz et alJ 
12001 " with the asymptotic HI velocities nearest our point- 
ings at ~510 km/s and ~'210 km/s respectively. We ob- 
served under topocentric radial velocities of -12 km/s and 3 
km/s toward UGC 7321 and UGC 1281 respectively. There- 
fore, we expect unresolved Ha emission at 6573.7±0.4A and 
6567.5±0.8A using the asymptotic values just quoted, in the 
observed frames of UGC 7321 and UGC 1281 respectively. 
The gas velocity dispersions in the 21 cm data are of the order 
7 km s^^. The 21 cm rotation curves change by ±10 km s^^ 
over our fields. These two values form the expected wave- 
length range, and the flux limits around these lines are flat to 
±100kms-i. 

Background galaxies produce the dominant flux in a num- 
ber of fibers. This is evident where we can measure red- 
shifts through emission lines identifiable as either Lyman- 
a, [OII]A3727, H/?, [OIII]A4959, or [OIII]A5007. For most 
of the background systems with emission lines the redshift 
is determined by the pattern of multiple emission lines. If 
the background galaxies have smooth continuum through our 
wavelength of interest their removal is accomplished in the 
continuum removal step. However, the possibility of spectral 
structure in the continuum across the corresponding Ha wave- 
length range leads us to mask those regions. Operationally, 
we mask a fiber if it displays a 5a significant value in its con- 
tinuum as estimated across all available wavelengths under 
inverse variance weighting. It is also possible that weak con- 
tinuum is coming from objects in the halo of the target galax- 
ies, in which case the desirability of a mask is less certain. We 
have performed all the emission line searches and limits with 
and without this masking process and found no detections in 
either case. The values we present as limits were made with 
the masks applied. 



3.5. Data co-addition and limits 

We show the derived limits in Table |2] We find no sig- 
nificant emission lines within the vicinity of the galaxies' 
velocities (defined as within ±500 km s^^) in any individ- 
ual fiber We next mask out continuum sources and apply a 
circular spatial filter as a 2D Gaussian function kernel with 
FWHM=10". Again, we find no significant emission. Fi- 
nally, we stack all fibers for which the model of i l2. 11 predicts 
/.J > 10^^^ erg/s/cm^/n". The choice of the cut in /i is not rig- 
orously determined, but judged as sensible from the shape of 
the surface brightness distribution in Figure[3]which indicates 
that many dozens of fibers contain predicted surface bright- 
nesses at roughly one-half the peak value. The co-addition 
of these fibers to one peak fiber will obviously yield an im- 
proved S/N. Different choices in the cut will lead to slightly 
different formal limits, but the fractional effect is small once 
large stack sizes of several hundred fibers are reached. The 
models used to select those fibers are those presented in Table 
[T]with an assumed F = 4 x 10~^^ s~^ and /? = 1.8. We use 
the nearly (Figure [U linear scaling between F and /i to deter- 
mine the true value of F. The models predict such averages to 
yield 1.7x10-1^ erg/s/cm^/D" for UGC 732 1 and 1 .8 x lO" 
erg/s/cm^/n" for UGC 1281. We again find no significant 
emission in the stacked spectra. These emission line searches 
were performed solely with errors base d on Poisson noise 
statistics and yielded no detections. In 0.61 we discuss ad- 
ditional systematic errors that degrade the final limits derived 
from purely Poisson errors in the data. By the models, the 
peak Ha surface brightness would have fallen in our fields 
for UVB strengths from 2xl0"i'* s^^ < F < 2x10-^2 g-i 
and warps of < 12.4° in UGC 7321 and4xl0-i^ s^^ < F < 
3x10-1'"^ s^i and warps of < 15.8° in UGC 1281. However, 
a radial displacement of the field would still give significant 
flux as seen in the contour plots, so we do not expect mis- 
alignments to affect the final limits. Figure [3] shows the posi- 
tions of the observations relative to several key features. The 
21 cm data contours are overlayed, the locations of masked 
background galaxies are shown, and the expected spatial pro- 
files of Ha emission is shown. We show in Figure H the sky 
spectra and the three types of spectral stacks to background 
subtracted data in UGC 7321. In Figure |5] we show the corre- 
sponding ones for UGC 1281. In neither case do we make a 
significant detection in Ha. 

Our selection of fibers for co-addition based on an assumed 
value of r leads us in turn to a lower limit on F. This may 
in principle introduce an error into our determination of F. 
However, both a rough estimation and then a detailed anal- 
ysis show that the non-linearity in this operation is negligi- 
ble. First, one measure of the spatial scale of the Ha surface 
brightness profile is the threshold radius, Tc, of Equation |5] 
Since Vc scales only as the natural logarithm of F, there is 
little change over the range of possible UVB strengths that 
we consider. The shape of the Ha surface brightness pro- 
file is also broad and smooth, from Figure [3] relative to the 
possible range of Vc- By selecting wide swaths of fibers for 
co-addition, the problem is particularly well behaved. Sec- 
ond, we verify these arguments with a numerical example. We 
simulated the surface brightness profiles for F = 8 x 10~^^ 
s^^, or five times lower than the nominal modeled value. 
We sampled the same set of fibers for co-addition as with 
the previous analysis. The model, average surface brightness 
was fi = 3.1 X 10^20 erg/s/cm^/D" and // = 3.9 x lO^^o 
erg/s/cm^/D" for UGC 7321 and UGC 1281 respectively, or 
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only 6% lower and 10% higher than the linear prediction. We 
conclude that the selection of co-added fibers based on the 
nominal UVB strength has negligible impact our final limit. 

3.6. Error assessment 

There are several potential sources of systematic error to 
the presented spectra. We have already discussed the uncer- 
tainties in the model-bas ed co nversion of Ha surface bright- 
ness to UVB strength in i l2.2l The uncertainty in the absolute 
spectral flux calibration due to the applied atmospheric extinc- 
tion curve is discussed in f|3] The uncertainty in the absolute 
spectral flux calibration due to the standard star observations 
is discussed in f|TT] We now analyze a final systematic re- 
garding the relative error determinations in the Ha spectra. 
We observe that the propagation of the errors from the data's 
original read noise and shot noise does not fully account for 
the variation in sky subtracted spectra. This is especially true 
under bright skylines. We discuss three possible causes with 
a focus on the variation of spectral resolution across different 
fibers. In any of the cases, the form of the systematic error will 
be to add a small percentage of the continuum subtracted sky 
background spectrum applied linearly with the random error. 

First, the instrumental spectral resolution varies by at most 
5% in different fibers due to small but detectable optical dis- 
tortions in the camera. We further measure from arc lamp 
exposures that the variation is 2.5% between the sky and sci- 
ence fibers in the UGC 7321 data and 1.5% in the UGC 1281 
data. These factors are presented in column 3 of Table |2] and 
scaled by the background subtracted sky spectrum and ap- 
plied as systematic errors in the spectra presented in Figures 
|4]and|5] This form of the systematic, as the fractional error in 
the dispersion times the background subtracted sky spectrum, 
can be derived simply by taking the first order expansion of 
a Gaussian function near the line center. Second, the fiber- 
to-fiber throughput can vary slightly between flat field cali- 
brations. The relative fiber-to-fiber throughput is calibrated 
with sky flats taken at dawn and dusk. This relative through- 
put has been measured to be stable to <5% over most nights. 
However, we find a maximum 15% fiber-to-fiber throughput 
variation in the UGC 7321 data due to poor fiber cable coiling 
practices. This error is very evident in the broadband estimate 
per fiber as shown in Figure |3] The error is less important for 
a continuum subtracted spectral element where most of the 
fiber-to-fiber throughput error subtracts out. The form of the 
throughput variation is that a few fibers experience a change 
with time, but the majority stay stable. We measure the rms 
throughput variation between all fibers to be far below 1%. 
Third, sky lines may vary across the ^ 1' separating the sky 
and science fibers. The UGC 7321 fiducial "signal" is well 
resolved from all known sky lines and only near OH lines, 
but the UGC 1281 "signal" is unresolved from an OH line 
and near the geocoronal Ha. Variations on such small spa- 
tial scales have not been observed, and the data are averaged 
over very long integration times and large ranges in zenith dis- 
tance, so we do not expect sky variation over our field-of-view 
to be a dominant error term. It is possible that the geocoronal 
Ha emission may vary within ~ 1', but small-scale variation 
is less likely for OH. We choose to parameterize the total ef- 
fect of these systematics in a conversative manner by deriving 
from the data themselves the systematic error based on the 
measured levels of spectral resolution variation seen between 
fibers. 

This systematic error strongly affects the UGC 1281 data 
since the lower resolution allows blending of night sky lines 



at the expected wavelength of Ha, but it is a less important 
component to the UGC 7321 error budget. As data from more 
fibers are coadded, this systematic error takes on greater im- 
portance in relation to the random error. We assess the 
distributions across 6300-6600A in each co-addition case in 
Table |2] The are simply calculated against a flat, zero flux 
line and can be visually judged in Figures |4] and |5] The dis- 
tributions look very symmetric around zero, and the reduced 
values are consistent with noise. The proper values val- 
idate our systematic noise estimates empirically. In fact, the 
additional noise estimates may be slightly conservative. One 
can visually note from Figures |4] and |5] that the values are 
even lower than the degrees of freedom {v) in the most impor- 
tant regions near the target wavelengths. 

3.7. Internal galactic extinctions 

Internal extinctions in disk galaxies at these scale lengths 
are very uncertain despite being a subject of detailed re- 
search (e.g. Byunetal. 1994). Matthews et al. (1999) see 
in UGC 7321 an abrupt truncation of resolved dust clumps 
beyond r«80" and fit a model of radially declining dust 
where, for our position around 250", there is no internal 
extinction. We have taken short VIRUS-P exposures cov- 
ering Ha and H/3 on the galaxy centers to derive conver- 
sative internal extinction upper limit estimates after correc- 
tion for Galactic extinction. We did not take deep enough 
exposures to measure accurate stellar populations and pho- 
tospheric Balmer absorptions ourselves, so we have re- 
lied on literature values appropriate to late type galaxies. 
From the Balmer decrements we measure A^fQ=-0.03±0.09 
magnitudes for UGC 7321 and A^fQ=-0.02±0.1 1 magni- 
tudes for UGC 1281 under the assumption that the absorp- 
tion equivalen t widt hs satisfy EW(Ha)n i,..;=EW(H/j)„ft,s=2A 
(iMcCaU et al.lll985l: ICalzetti et all 1 1 994 . As the extinction 
estimates are consistent with zero, we apply no dust correc- 
tion to our results. 

4. DISCUSSION 

The flux decrement method is currently the most widely 
used method to estimate the UVB strength at high red- 
shift. Under the fluctuating Gunn-Peterson approximation 
dCroft et al.llI998l) . the Lyman-a forest optical depth distri- 
bution should have a normalization that depends only on well 
constrainted cosmological parameters and the UVB strength. 
The IGM temperature and density distributions may have 
some systematic uncertainties that propagate into knowledge 
of the UVB, but they are not likely the leading uncertain- 
ties. The more likely dominant uncertainties in flux decre- 
ment modeling are the source emissivities. At z ^ 1, the 
Lyman limit mean free path becomes larger than the horizon, 
so the UVB strength at z=0 is influenced by source evolution 
across this redshift range. AGN and stellar population lumi- 
nosity functions, both observed and modeled, generally agree 
to better than an order of magnitude over these redshifts. The 
least constrained input to flux decrement modeling is the es- 
cape fraction for ionizing photons in galaxies, particularly at 
low redshift and low luminosity. We believe our measurement 
is best interpreted as an indicator of a low escape fraction. 

Our most constraining (5(t) spectral limits are F < 1.7 x 
10-14 s-i in UGC 7321 and F < 13.5 x IQ"^-^ s'^ in 
UGC 1281 again assuming /? = 1.8. Several benchmarks, 
both empirical and theoretical, exist with which to compare 
these limits. Figure|6]shows the UVB strength against redshift 
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determined by many gr oups. The lowes t redshift proximity 
effect limit comes fro mlKulkarni & Falll (Il993h with analy- 
sis of 13 quasars from lBahcall et alj ( |1993^ between 0.16 < 
z < 1.00 at r(z = 0.5) = 2.0tl% x 10"" 8"^ However, 
the proximity effect method h as been shown to have a high 
bias that depends on halo mass dPaucher-Giguere et al.l2008a) 
an d should be interpreted with c are. The theoretical model 
of iFaucher-Giguere et al.l (l2009h gives a drop in the UVB 
strength by a factor of 3.4 between z=0.5 and z=0.0 leav- 
ing this measurement consistent with our current limit. This 
agreement is interesting and somewhat unexpected given the 
bias of proximity effect measurements. The only existing low- 
z flux decrement hmi t is T{z ^ 0.17) = S.Oifo x 10"^'' s"! 
(iDave & Tripp! |200]T) . The theoretical model itself, normal- 
ized by the flux decrement method, predicts r(z = 0) = 3.8 X 
10"^"' s^^ which is much higher than our new limit. There ex- 
ists a second set of unpublished theor etical predictions from 
F. Haa rdt and P. Madau discussed in iFaucher-Giguere et al.l 
(|2009) giving r(z = 0) = 1 x IQ-^^ ^-i -j-jjg j^jfter model 
used a constant 10% escape fraction of ionizing photons and 
an unspecified star formation history while the former used 
a compl etely theoretical and simula tion-based star formation 
history (iHernquist & Springelll2003T) and a scaling of the stel- 
lar UV emissivity based on high redshift flux decrement mea- 
surements that contains the escape fraction. A comparison 
to Lyman-break galaxy (LBG) lumin osity functions led that 
gro up to require only fesc.abs ~ 0.5% (IFaucher-Giguere et al.l 
l2008bi) . The direct measurement of galactic escape frac- 
tions is difficult due to the low values involved. While UV 
bright samples c an range up to » 3% in absolute Lyman limit 
escape fraction dShaplev et al.l |2006|). a presuma bly lower- 
mass sample yielded (2 ± 2)% dChen et al .1120071) . Theoret- 
ical work shows a stron g decrease i n f^sc with star forma- 
tion rate and halo mass (iGnedin eT al. 2008) below Mtot ~ 
IO^^A/q, and lower redshift observations of popu lations sim- 
ilar to LBGs show a potential redshift evolution (ISiana et al.l 
I20T0I) with fescMhs < 0.8%. There IS no reason yet to sup- 
pose a lower bound to the escape fraction. If we interpret 
our li mit as a scaling of the escape fraction from the mod- 
els in IFaucher-Giguere et al.l ( l2008bl) at low redshift, we find 

fesc,abs < 0.2%. 

It is unlikely that systematics from the model assumptions 
in our analysis can cause the disagreement. Contaminating 
ionization from the galaxies' forming stars would bias our 
measurement high, only making the disagreement more se- 
vere. We further note that the degree of contamination can be 
measured by anomalous [NII]A6548 to Ha ratios (BFQ) and 
should not, in principle, limit this type of measurement. There 
has been a large body of work on low strength star formation 
beyond the optical radii in local galaxy disks, usually labelled 
extended UV disks (XUV), fostered by far UV (FUV,1350- 
1750A) and near UV (17 50-2750A) Galaxy Evolution Ex- 
plorer (GALEX) data (e.g. lThilker et al.ll2007l) . Narrowband 
Ha imaging and spectroscopy have revealed that ^ 10% of 
gas rich disks (Werket al. 2010a b; Herbert-Fort et al. 20l3) 
host outlying Ha emitting complexes as either compact HII 
regions or dwarf satellite companions. The common Ha 
fluxes observed so far are of the order of a few times 10^^^ 
erg/s/cm^. Any such systems would have been found in our 



data as strong detections limited in size to a few fibers. The 
expectation of large-scale, diffuse UVB Ha emission should 
discriminate reliably against compact XUV Ha emission. We 
have also visually inspected the target galaxies' GALEX data 
which have not yet been analyzed in any XUV focused work. 
UGC 1281 has only been covered in the rather shallow all- 
sky survey mode. UGC 7321 has been covered for 2.8ks in 
the NUV and 1 .7 ks in the FUV under guest investigator cy- 
cle 4 pr oposal ID 095 (PI: J. Lee) as part of the IIHUGS 
project (iLee et al.ll2009h . In neither system is there evidence 
for an extended UV disk beyond the DSS2-red^ limiting con- 
tours. Finally, these contamination issues are speculative un- 
til a putative UVB Ha detection is made. The only possible 
systematics that could have made a low bias to our limit are 
unaccounted for dust or gas distribution parameters, such as 
inclination, far beyond the range we have constrained. 

We have made our first analysis under the assumption that 
the gas distribution extends beyond the HI data limits with 
the same exponential form as at smaller radii. This assump- 
tion, motivated by the thin and regular HI distributions and 
lack of nearby companions, has the strongest impact on our 
interpretation. An alternative estimate without this assump- 
tion, taking only fibers that overlap with the observed HI sig- 
nal, yields a very comparable limit of F < 2.3 x lO"^"' s~^ 
at 5(7 significance in UGC 7321. This agreement essentially 
comes about because our original model predicts only a minor 
Ha contribution at the discarded positions under the modeled 
UVB strength. Nevertheless, there is no reason to assume 
the presence of an HI edge since the radio observations de- 
tect the gas up to the column densities where the sensitivity 
runs out. This result raises the question whether a redshift- 
dependent escape fraction is manifesting in galaxies. Alterna- 
tively, our new limits may be saying that the UVB strength, 
as estimated through flux decrement measurements, has been 
overestimated at all redshifts. The latter choice would upset 
the apparent agreement between current models and reioniza- 
tion constraints. Either case will require some modification 
to the UVB strength model and its implementation in struc- 
ture formation simulations. We intend to pursue our measure- 
ments of these and other superthin galaxies to greater depth in 
order to arrive at a detection of r(z = 0). 
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APPENDIX 

FULL SOLUTION TO THE GENERAL Ha SURFACE BRIGHTNESS 

We give here the derivation of the general Ha surface brightness at field positions 61 and 62. The special case for 6i=rc and 
62=0 was derived as Equation |9] That case is simplified since the line of sight integration can proceed from zero to infinity 
without intersecting the photoionization boundary and has symmetry between positive and negative values of p. For the general 
case, the simple task presented in this Appendix is to find the possible geometrical intersections of z from Equation [T] and Zc{R) 
from Equation |6] as a function of p under inputs i, bi, and 62- This may have zero or two intersections labelled as pri and pr2- 
Once found, the general expression for /x then follows Equation lAll 

fJ'[t>l,D2) - ^ 



jP.^ exp(_ 2PCOS. _ h/^^^)dp+ (Al) 
I /£exp(-^ - : roots in p 



with 7 



The first necessary condition for any intersection to occur is evidently expressed in Equation lA2l as the largest possible distance 
for an intersection to lie from the galaxy center is Vc while the closest possible approach for a sight line is bi. 

hi < Tc (A2) 

The intersections in p can be expanded into simple quadratic equations. Each of the two potential roots from the quadratic 
solution is double valued when considering intersections with both signs of the Zc{R) surface leading to four possible roots. 
However, only at most two of the roots will be physical with the rejected two lying on extrapolations of Zc{R) at R{p) > Vc or 
\z{p)\ > Tc X hz/hr- The intersections with these surfaces lead to possible limits psi, Ps2, Pss, and ps4 expressed in Equations 
IA3HA6I The most constraining limits are then the values between these four with the smallest absolute values leading to Equations 
IA7HA8l for the active limits pn and pi2. 

-62 sin I - 

Psi = (A3) 

cost 



Ps2 = '^^-^^-V^l ^A4^ 



-62 sini + jfvc 

Ps3 = (A5) 

cost 



b2 cosi+ ^rl - b\ 

Psi = r^. (A6) 

pii = ma.x{psi,ps2) (A7) 

P12 = "cainipss, Psa) (A8) 
The potential intersections with Zc{R) can be directly evaluated as ppi, pp2, Pp3, and pp4 as given in Equations IA9II A 1 2l 



(tj^) ''c '^'^^ * ^ (f^) ^2 sill i COS 1 — 62 sill i cos i + 
Pp^ ^ 7~T2 



cos^ i — sin^ i 



/3pi = 2 I — I 60 cos^ i sin^ i — 2 I — ) 62''c cos^ i sini + 5? [ — I cos^ i 
\hz J \hz J \hz J 

+bl ' sin2 t + rl sin2 t + ' bl sin^ z - 2 sin^ i - b\ sin^ i (A9) 

(r^) * ~ (x^) ^2 sin i cos z — ^2 sin i cos i — ^ Ppi 



Pp2 



■ 2 ■ 

sm I 



3pi = 2 ( ^ ) b\ cos^ i sin^ i - 2 ^s^-c cos^ i sinz + h{ '^'^^ * 

+&1 ( jl^^j Sin2i + r2sin2i+ ^j^^j blwi" i - 2{^^r^b2^u^ i - b\€Ya^ i (AlO) 
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Jp2 



Pp3 



(ti^) '^'^^ * ^ (r^) ^2 sin i cos i — 1)2 sin i cos i + 



bi cos" i sin" i + 2 



cos^ I — sm I 



biTr cos^ i sin i + b\ . 



cos^ i 



sin^ i + r'i sin^ i 



&2 sin"* i + 2 ( ^ ) rc62 sin^ i - bf sin^ « 



(All) 



(ti^) '■^'^^ * ^ (ti^) ^2 sin i cos i — ^2 sin i cos i 
^p4 = ; r2 



^ I cos2 i — sin^ i 



3p2 = 2 ( ) ^2 cos^ i sin^ i + 2 ( 1 62''c cos'^ i sin z + 5^ ( ^ ) cos^ i 



+bl (^'^^ ' sin^ t + rl sin^ z + (j^^ ' 6^ sin^ i + 2 r,62 sin^ i - bf sin^ i (A12) 

The comparisons to the limits pn and p;2 discard unphysical values in Equations I A 1 3 ||AT41 where the final limits of integration 
are found. 

Pri = min(x e {ppi,Pp2,Pp3,Pp4:} ■ Pii < X < P12}) (A13) 
Pr2 = max(x e {/0pi,/3p2,/3p3,/0p4} : pn < x < P12}) (A14) 
With the integration boundaries now well defined, p{bi, 62) can easily be obtained through numerical integration. 
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Figure 1. The fractional eiTor in a linear relation between the Ha peak surface brightness and the UVB photoionization rate, F, under different UVB strengths. 
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Figure 2. The single position angle, parameterized fits to the HI distributions shown as solid lines. Because we assumed a single position angle and single radial 
scale length, the ranges over which we fit the HI distributions must be somewhat controlled and limited to the radii near the Ha observations. The predicted Ha 
surface brightness profiles are also shown against the right side axis as dotted lines. The horizontal arrows indicate the correct axis for each profile. The breaks 
at large radii in the 21 cm profiles are due to intersections with the UVB photoionization fronts. Top The data along the major axes. Bottom The data along the 
minor axes at particular offsets. Left Cuts along the midplane and normal to it offset by 165" with data, models, and 68% confidence intervals in UGC 7321. We 
restrict the fit to points > 160" from the galaxy's center as indicated by the vertical dotted line. Using the best parameters from Table [T] the threshold radius 
(Equation|5) with the nominal value of F = 4 X lO^^^ s^^ is predicted to be at = 13.4kpc. Our spectroscopic data cover regions from R= 9.5kpc to 
R= 14.6kpc. Right Similarly, data and fits to UGC 1281. The offset here is 200" from the minor axis. The points between the dotted lines at 160" and 220" 
form the restricted range of the fit as a 8° warp becomes important beyond. This fit appears poorer because of the larger warp, but a fit to allpoints at R > 220" 
returns the same radial scale length to within the Monte Carlo eiTors. Using the best parameters from Table[T] the threshold radius (Equation[5) is predicted to be 
at rc = 7.4kpc. Our spectroscopic data cover regions from R= 5.7kpc to R= 9.1kpc. Since we only show one dimensional cuts, these figures do not show all 
the datapoints used in the fits. 



Table 1 

HI based model parameters and Ha surface brightness predictions* 



Galaxy 


no 
(cm-3) 




hz 
(pc) 


h, 

(kpc) 


i 

(°) 


PA 

(°) 


A»0 

t 


t 


tt 


mi 


UGC 7321 
UGC 1281 


T -1+3.5 
-'■-'-1.7 


426. 
303. 


+ 120. , 
-88. "10 
+ 70. J 
-58."5 


2.i2tE;:?irfio 

, n+0.19 J 


82.810-1 
84.9tt;3 


-lOO.litO.l 
-141.3±0.3 


18.4 
21.4 


l°-'-0.7 

l9Atli 


n.9tl-l 


22.5it;l 



* Fit under the restricted radial ranges shown in Figure |2]assuming r = 4 X lO^i-* s^^ and /3 = 1.8. 
t 10"^" erg/s/cm^/n" 

-f 10~^° erg/s/cm^/n", .smoothed by a circular 10" FWHM kernel 

10"^° erg/s/cm^/n' ', average for all fiber positions with predicted values of fi >10 erg/s/cm^/n" 
** lO"^'^'^ erg/s/cm^/n", based on the HI bounded area (Equations 6-8 ot lWevmann et alJ|20Ql|) ) 
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Table 2 

Error budget and limits to the UVB strength 



Co-addition 


Poisson 


Resolution 


Flux 


SB 


Model 


r{z = 0) 




type 


error 


systematic 


calibration 


upper limit 


systematic 


upper limit 


(1) 


(2) 


(3) 


(4,%) 


(5) 


(6,%) 


(7) 


(8) 



UGC 7321 each single fiber 


11.0 


0.8 


8.9 


64 


+5.4/-4.9 


15 


353/468 


UGC 7321 smoothed 


2.8 


0.4 


8.9 


17 


+6.6/-4.2 


4.4 


257/454 


UGC 7321 radio bound stack 


1.8 


0.4 


8.9 


12 


+19.3/-7.9 


2.3 


545/462 


UGC 7321 full stack 


0.9 


0.4 


8.9 


7.1 


+6.0/-1.8 


1.7 


497/454 


UGC 1281 each single fiber 


18.6 


29.5 


4.3 


250 


-I-57/-13 


53 


68/136 


UGC 1281 smoothed 


6.6 


8.9 


4.3 


81 


+29/- 12 


19 


218/136 


UGC 1281 radio bound stack 


6.0 


8.9 


4.3 


78 


-1-46/- 15 


27 


261/134 


UGC 1281 full stack 


2.0 


8.9 


4.3 


57 


-H9.5/-6.1 


14 


50/134 



Detection and model method. Smoothed refers to a 10" FWHM Gaussian, circular kernel. 
Icr (lO""^*^ erg/s/cm^/n") in the spectral data from Poisson noise. 

Icr (10"^*^' erg/s/cm"^/n") in the spectral data from spectral resolution or sky line variation. See ^3.6| 
1 (T (%) flux calibration systematic including 2. 1 % for airmass/extinction error. 

5(7 (lO"^'^ erg/s/cm^/n") limit in surface brightness. No detections of significance were found, so this limit results simply from summing 
columns 2 and 3, multiply by one plus the percentage in column 4, and finally multiplying by five. 

1(T (%) model surface brightness systematic. These values are derived from the Monte Carlo tests of ^2.2| 
'"^^ 5a- (10~^^ s~^) total limit assuming /3 = 1.8. The achieved Ha surface brightness limit is compared to the low bound of the modeled 
Hq surface brightness to create this final, linearized estimate from the modeled value ofF — 4 X lO^^'^ s^^. This limit results simply by 
multiplying column 5 by one plus the percentage in column 6, multiplying by the t 
lower bound to either column 7, 8. 9, or 10 in TablefTldepending on the limit type. 

of all pixels between 6300-6600A in the spectrum and the degrees of freedom. 
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Figure 3. Left Reconstructed VIRUS-P continuum image of the UGC 7321 outskirts centered at qj2ooo=12:17:16.4 and 5j2000=+22:31:33 or sd250" off the 
minor axis. The continuum estimation is made through the entire available spectral range from 6040-6740A with the colorbar in units of 10~^^ erg s~^ cm~^ 
A~^. The dark, circled objects are masked as background galaxies, many known to be background by their emission lines at redshifts high er th an the target 
galaxy's redshift. One can see some broad structure in the continuum map due to small residuals in the fiber-to-fiber throughput as described in i|3.6l especially in 
the UGC 7321 data. The green contours trace the HI column densities in steps of (10,19,36,67,126,238,448,845) xlO^* cm^-^. The red, more extended contours 
trace the predicted Ho surface brightness assuming r = 4 X 10" s"^ and (3 = 1.8 in contour levels of (0.1,0.24,0.57,1.4,3.3,7.9,19) X 10"^" erg/s/cm^/D". 
The two innermost red contours enclose the surface brightness maxima. Positions closer to the center again become dimmer in Ha since portions of the gas, in 
projection, stay neutral at smaller radii. The fibers used in sky subtraction are all those outside the second outermost red contour. We draw the second, seventh, 
and eighth contours thickly to highlight these regions. As a scale reference, the fiber diameter is 4'' 1. Right The same display for UGC 1281 with central position 
c*j2000=l:49:15.8 and (5j2ooo =+32:3 1:46 or 300" off the minor axis. The continuum estimation is made through the entire available spectral range from 
4700-6990A. Here, many more background galaxies are found. In UGC 1281, we took data at two overlapping fields. The central positions covered by both 
square pointings have the best depth. 
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Figure 4. Selected spectra around Ha in UGC 7321 presented in units of surface brightness. The expected wavelength for emission is shown with a dotted line. 
The frames from top to bottom show the background sky, the background subtracted spectrum for a typical fiber that does not display continuum, the spectrum at 
the same position after being smoothed by a 10" FWHM circular Gaussian kernel, the data bounded by HI signal, and finally the stack of the 358 fibers predicted 
to be the brightest by the model. The errorbars consist of the Po isson , observational error and the systematic spectral resolution error of columns two and three 
in Table|2]only. The spectral resolution systematic, discussed in i|3.6l is most important under the bright skylines and does not dominate at the target wavelength. 
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Figure 5. Selected surface brightness spectra around Ha in UGC 1281. The format is the same as in Figure|4] In this case, 313 of the brightest expected fibers 
form the final stack. 
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Figure 6. A compilation of photoionization rates across redshift. Most of the literature compilation s come from Table 2 in 'Faucher-Giguere et al.' ( 2008^ 
and Table 1 in Faucher-Gigue re et al.l )2008bf) . The flux decrement measurement at z ~ 0.17 is from lDave & TrippI (200 10. The low redshift, H« limit from 
[Wevmann et al. ( 2001 ) (2a) has been the deepest z=0 limit before this work. The UVB fitting function comes from|Fardal et al. ( 1998) and the newer simulation 
from Faucher-Giguere et al. 1 200^. Our work's new limit is well below the flux decrement normalized simulation and challenges one or more of the model 
assumptions. Some points have been slightly shifted in redshift for visual clarity. 



